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ABSTRACT 

Star formation activities in the NGC 2024 FIR 4 region were studied by imaging centimeter continuum sources 
and water maser sources using several archival data sets from the Very Large Array. The continuum source VLA 9 is 
elongated in the northwest-southeast direction, consistent with the FIR 4 bipolar outflow axis, and has a flat spectrum in 
the 6 .2-3.6 cm interval. The three water maser spots associated with FIR 4 are also distributed along the outflow axis. 
One of the spots is located close to VLA 9, and another one is close to an X-ray source. Examinations of the positions of 
compact objects in this region suggest that the LIR 4 cloud core contains a single low-mass protostar. VLA 9 is the best 
indicator of the protostellar position. VLA 9 may be a radio thermal jet driven by this protostar, and it is unlikely that 
LIR 4 contains a high-mass young stellar object (YSO). A methanol 6.7 GHz maser source is located close to VLA 9, 
at a distance of about 100 AU. The LIR 4 protostar must be responsible for the methanol maser action, which suggests 
that methanol class II masers are not necessarily excited by high-mass YSOs. Also discussed are properties of other 
centimeter continuum sources in the held of view and the water masers associated with LIR 6 n. Some of the continuum 
sources are radio thermal jets, and some are magnetically active young stars. 

Subject headings: ISM; individual objects (NGC 2024) - ISM: jets and outflows - ISM: structure - masers - stars; 
formation 


1. INTRODUCTION 

The NGC 2024 molecular ridge is a filamentary cloud ex¬ 
hibiting various star formation activities (Mezger et al. 1988; 
Visser et al. 1998; Matthews et al. 2002), at a distance of 415 
pc from the Sun (Anthony-Twarog 1982). The molecular ridge 
consists of dense cores LIR 1-7, and they contain the youngest 
objects in the NGC 2024 region (Mezger et al. 1992; Wilson 
et al. 1995; Chandler & Carlstrom 1996). Mezger et al. (1988, 
1992) suggested that the LIR 1-7 cores are isothermal protostars 
without luminous stellar objects, but later studies showed that 
some of them are associated with young stellar objects (YSOs) 
and molecular outflows (Chandler & Carlstrom 1996). While 
the structures and star formation activities of LIR 5 and LIR 
6 are known relatively well (Richer 1990; Richer et al. 1992; 
Wiesemeyer et al. 1997; Lai et al. 2002; Alves et al. 2011; Choi 
et al. 2012b), the nature of the other LIR cores are not clearly 
known. Recent detections of molecular masers and X-ray emis¬ 
sion suggest that LIR 4 is yet another active site of star forma¬ 
tion, showing interesting physical phenomena at small scales 
(Minier et al. 2003; Skinner et al. 2003; Choi et al. 2012a; Green 
et al. 2012 ). 

The dense core LIR 4 contains a low-mass protostar and 
exhibits star formation activities such as a reflection nebulos¬ 
ity and a molecular outflow (Moore & Chandler 1989; Lis et 
al. 1991; Moore & Yamashita 1995; Chandler & Carlstrom 
1996). By contrast, Minier et al. (2003) detected a CH 3 OH 
class II maser source and suggested that LIR 4 contains an 
intermediate/high-mass YSO. It is unclear whether or not the 
outflow-driving protostar and the CHsOH-exciting YSO are the 
same object. High-resolution imaging is needed to identify the 
YSO in the LIR 4 region and to investigate its nature. Interfer¬ 
ometric observations in the radio continuum can increase our 
understanding of the LIR 4 region, but the imaging of compact 
sources is hampered by the strong emission from the extended 
H II region (Crutcher et al. 1986; Gaume et al. 1992). The ef¬ 
fects of the extended emission can be reduced by excluding the 


visibility data of short baselines (Rodriguez et al. 2003; Choi et 
al. 2012 b), but the resulting images contain artifacts introduced 
by the Altering and need to be analyzed carefully. 

The 1.3 cm H 2 O maser of low-mass protostars is a rarely de¬ 
tectable phenomenon (Kang et al. 2013). Once detected, how¬ 
ever, the maser indicates the existence of shocked gas and can 
provide important information about the star formation activi¬ 
ties of YSOs driving the shock (Genzel & Downes 1977; Elitzur 
1992; Luruya et al. 2003). Particularly, the small source size and 
high intensity allow interferometric observations, which is es¬ 
sential in the investigation of small-scale structures. The H 2 O 
maser associated with LIR 4 was first reported by Choi et al. 
(2012a). They used a single-dish telescope, however, and the 
large beam size made it difficult to understand the relation be¬ 
tween the H 2 O maser source and other compact objects in the 
region. 

In this paper, we present the results of archival observations 
of the NGC 2024 LIR 4 region in the centimeter continuum or 
in the H 2 O maser line with the Very Large Array (VLA) of the 
National Radio Astronomy Observatory (NRAO). The archival 
data sets are described in Section 2. In Section 3, we report the 
results of the radio imaging. In Section 4, we discuss the star¬ 
forming activities of LIR 4. In Section 5, we describe other ob¬ 
jects in the held of view. A summary is given in Section 6 . 

2. DATA 

The VLA data sets presented in this paper were retrieved 
from the NRAO Data Archive System in the form of raw tele¬ 
scope data files. Table 1 gives a summary of the data sets. The 
calibration and imaging were done using the Astronomical Im¬ 
age Processing System software package of NRAO. 

2.7. Centimeter Continuum 

We searched VLA data sets in the NRAO Data Archive Sys¬ 
tem for continuum observations of the LIR 4 region to investi¬ 
gate the nature and structure of the continuum source associated 


1 


Table 1 

Summary of the VLA Data Analyzed 


Project 

Observation Date 

Configuration 

Data Type 

Fph'' 

AG 421*' 

1994 Aug 13 

B 

2.0, 3.6, and 6.2 cm continuum 

0.97 Jy at 3.6 cm 


1994 Oct 24 

C 

2.0 and 3.6 cm continuum 

1.00 Jy at 3.6 cm 


1995 Aug 17 

A 

2.0, 3.6, and 6.2 cm continuum 

1.30 Jy at 3.6 cm 

AF354 

1999 Feb 26 

DnC 

1.3 cm H 2 O line 

~1.73 Jy at 1.3 cm 

AR465 

2002 Mar 2, 3, 8 

A 

3.6 cm continuum 

0.76 Jy at 3.6 cm 

AM 749 

2003 Mar 14 

D 

1.3 cm H 2 O line 


AM 780 

2003 Sep 15 

A 

1.3 cm H 2 O line 

0.83 Jy at 1.3 cm 


Flux density of the phase calibrator (QSO B0539-057), showing the time variability. 

The AG 421 data set presented in this paper is a subset of the full AG 421 data set (see Section 2.1.1). 


with the FIR 4 protostar. Two data sets were found to be useful. 

2.1.1. Project AG 421 

The NRAO observing project AG 421 includes data from 
observations made in various array configurations. Since we are 
mainly interested in compact sources, D-array and DnC-array 
data were ignored, and data from more extended configurations 
were considered. The B-array, C-array, and A-array observa¬ 
tions were made in 1994 August, 1994 October, and 1995 Au¬ 
gust, respectively. 

The NGC 2024 region was observed in the standard contin¬ 
uum modes of the C band (4.86 GHz or A = 6.2 cm), X band 
(8.44 GHz or A = 3.6 cm), U band (14.94 GHz or A = 2.0 cm), 
and K band (22.46 GHz or A = 1.3 cm). The phase-tracking 
center was aiggo = 05‘'39“14!30, Jiggo = -01°55'55"0, which 
corresponds to IRS 2. The phase was determined by observ¬ 
ing the nearby quasar 0539-057 (QSO B0539-057). The flux 
was calibrated by observing the quasar 0134H-329 (3C 48) or 
1328H-307 (3C 286). The flux densities of 0539-057, derived 
from a comparison of amplitudes, were 0.96/1.05/1.23 Jy at 6.2 
cm, 0.97/1.00/1.30 Jy at 3.6 cm, and 0.94/1.00/1.16 Jy at 2.0 cm 
for the B-array/C-array/A-array observations, respectively. The 
observations were made with relatively short scans. A typical 
on-source integration time in each array configuration was ~14 
minutes for each wavelength band. 

Maps were made using a CLEAN algorithm. To avoid the 
adverse effects of the strong extended emission associated with 
the nearby ionization fronts, the visibility data of short baselines 
(uv < 50 kA) were excluded. This filtering was applied to the 
data of all array configurations. For the 6.2 cm continuum, the 
C-array data were not used because all the baselines were short 
(uv < 53 kA). For the 2.0 cm continuum, a uv taper of 750 kA 
was applied to make the beam size comparable to that of the 3.6 
cm continuum. The 1.3 cm continuum data did not produce a 
useful map. With a robust weighting, the 6.2 cm, 3.6 cm, and 
2.0 cm continuum data produced synthesized beams of (y/51 x 
(y.'46, (y.'33 X (y.'29, and (y/33 x (XQO, respectively, in FWHM. 

As a result of the short integration time and the exclusion 
of short-spacing data, the images have relatively low qualities 
and contain numerous artifacts in a form of spurious intensity 
peaks. Detection of sources is limited by confusion with these 
unreal peaks. Therefore, we focused on measuring the flux den¬ 
sities of the sources at known positions (Rodriguez et al. 2003) 
and did not try to identify new sources. The AG 421 data set 
is useful for constructing the continuum spectra of relatively 
strong sources, and the AR 465 data set (see below) is useful 
for studying source positions and structures. 


2.1.2. Project AR 465 

The NRAO observing project AR 465 includes data from 
three tracks of A-array observations made in 2002 March. The 
NGC 2024 region was observed in the standard AT-band con¬ 
tinuum mode (8.46 GHz or A = 3.6 cm). The phase-tracking 
center was agooo = 05‘'41“44!90, Jgooo = -01°55'54"0. Details 
of the observations and results were presented by Rodriguez et 
al. (2003). The calibration was done in the way as described 
by Rodriguez et al. (2003). The total on-source integration time 
was 5.8 hr. 

The visibility data from the tree tracks were combined, and a 
map was made using a CLEAN algorithm. The visibility data 
of short baselines (uv < 50 kA) were excluded. With a ro¬ 
bust weighting, the 3.6 cm continuum data produced a synthe¬ 
sized beam of FWHM = O'/25 x (y.'23. (The map presented by 
Rodriguez et al. 2003 was made using data with uv > 100 kA, 
and the data from the three tracks were combined in the image 
space.) 

2.2. Water Maser 

We searched VLA data sets for HgO maser observations of 
the FIR 4 region to investigate the relation between the maser 
source and the FIR 4 protostar. Three data sets were found to 
have pointing centers close to FIR 4. In these observations, FIR 
4 is located near the edge of the field of view. 

2.2.1. Project AF 354 

The NRAO observing project AF 354 includes data from the 
DnC-array observations made in 1999 February. Details of the 
observations and results were presented by Furuya et al. (2003). 
For the HgO Gig 623 line (22.235077 GHz), the spectral 
window was set to have a velocity resolution of 0.33 km s“^. 
The phase-tracking center was aigso = 05^39™ 13!35, ( 5 ig 5 o = 
-01°57'19'.'0. The angular distance between the phase-tracking 
center and FIR 4 is ~70" (Figure 1). For comparison, the an¬ 
tenna primary beam has an FWHM = 122". The phase calibra¬ 
tor was the quasar 0539-057 (QSO B0539-057). The flux was 
calibrated by setting the flux density of 0539-057 to 1.73 Jy, 
which makes the flux scale consistent with that of Furuya et al. 
(2003). 

Maps were made using a CLEAN algorithm. With a natural 
weighting, the HgO line data produced a synthesized beam of 
FWHM = 3"1 X 2'.' 1. The imaging was done in a square box of 
160" X 160", which is larger than the size of the primary beam, 
in order to include FIR 4 in the imaging area. 
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Table 2 

NGC 2024 Continuum Source Parameters 


Source 

6.2 cm 

Flux Density'* 
3.6 cm 

2.0 cm 

o'* 

Associated 

Objects 

Source 

Type" 

VLA 1 

0.70 ± 0.04 

0.58 ± 0.03 

OFOV 

-0.33 ±0.17 


YSO (thermal jet) 

VLA 2 

0.47 ± 0.05 

0.17 ±0.02 

OFOV 

-1.90 ±0.33 


YSO (magnetic activity) 

VLA 3 

6.64 ± 0.05 

6.84 ± 0.03 

4.37 ± 0.04 

0.05 ± 0.08 


YSO (thermal jet) 

VLA 8 

2.04 ± 0.03 

1.30 ±0.02 

0.39 ± 0.02 

-0.82 ± 0.09 


Background object 

VLA 9 

0.51 ±0.03 

0.50 ± 0.02 

0.97 ± 0.02 

-0.04 ±0.14 

FIR 4 

YSO (thermal jet) 

VLA 16 

0.38 ± 0.03 

0.58 ± 0.02 


0.76 ±0.16 


YSO 

VLA 19 

13.29 ± 0.03 

24.82 ± 0.02 

39.92 ± 0.02 

1.13 ±0.08 

IRS 2 

Early-B type star 

VLA 21 

0.78 ± 0.03 

1.56 ±0.02 

OFOV 

1.26 ±0.11 


YSO (magnetic activity) 

VLA 24 

5.26 ± 0.03 

6.06 ± 0.02 

OFOV 

0.26 ± 0.08 


YSO (magnetic activity) 

VLA 25 

0.45 ± 0.04 

0.62 ± 0.04 


0.59 ± 0.21 

IRS 4 

YSO (thermal jet) 


Notes.—S ee Table 1 of Rodriguez et al. (2003) for source positions and a more extensive list of associated objects. 

^ Total flux densities in mjy, corrected for the primary beam response. The uncertainties are the statistical noise of each 
image. OFOV: located outside the field of view. 

^ Spectral index in the 6.2-3.6 cm interval: Fi, oc u°‘, where iT is the flux density, and u the frequency. The uncertain¬ 
ties include the statistical noise and the uncertainty of the flux scale (assumed to be 3%). 

Nature of the radio source, based on the spectral slope, size, variability, and infrared counterpart. See Section 5.1 for 
details. 
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Fig. 1.— Color composite image of the 3.6 (blue), 4.5 (green), and 
5.8 (red) /rm bands toward the NGC 2024 region from the Spitzer data 
archive. Plus signs: 870 ^m continuum sources FIR 1-7 (Mezger et al. 
1992). Source names are labeled on the right edge of the figure. Dashed 
circle: FWHM area of the VLA primary beam centered at the pointing 
center of the AF 354 observations. The brightest object in the infrared 
image is IRS 2 (VLA 19). 


2.2.2. Project AM 749 

The NRAO observing project AM 749 includes data from the 
D-array observations made in 2003 March. For the H 2 O line, 
the spectral window was set to have a velocity resolution of 1.3 
km s“^. The phase-tracking center was a 2000 = 05^41™44!765, 
1^2000 = -01°55'53'.'76. The angular distance between the phase¬ 
tracking center and FIR 4 is ^70". The phase calibrator was 
the quasar 0541-056 (QSO B0539-057), the same as that of 
AF 354. The flux was calibrated by setting the flux density of 
0541-056 to 0.828 Jy (see below). 

Maps were made using a CLEAN algorithm in the area as 
described above. With a natural weighting, the H 2 O line data 
produced a synthesized beam of FWHM = 4'.'3 x 3'.'4. 

2.2.3. Project AM 780 

The NRAO observing project AM 780 includes data from 
the A-array observations made in 2003 September. The velocity 
resolution, phase-tracking center, and phase calibrator were the 
same as those of AM 749. The flux was calibrated by setting 
the flux density of the quasar 0137H-331 (3C 48) to 1.281 Jy. A 
comparison of the amplitude gave a flux density of 0.828 Jy for 
0541-056. 

Maps were initially made in the area as described above, 
which revealed detectable emission near FIR 6n only. Final 
maps were made in a square box of 37" x 37" around the 
phase-tracking center. With a natural weighting, the H 2 O line 
data produced a synthesized beam of FWHM = O'.' 13 x O'.' 11. 

3. RESULTS 

3.1. Centimeter Continuum 

Several sources were detected in the centimeter continuum 
images of the AG 421 data set. Table 2 lists the flux densities 
of the sources detected in multiple wavelength bands. VLA 9 is 
the radio source associated with FIR 4, and Figure 2 shows the 
images in the three bands. The source position agrees with that 
of the 3.6 cm source reported by Rodriguez et al. (2003). 
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Fig. 2.— Maps of the 6.2, 3.6, and 2.0 cm continuum emission toward the FIR 4 region, from the AG 421 data set. The contour levels are 1, 2, 4, 
and 8 x 3a, where a is the rms noise of each image: 0.03, 0.015, and 0.018 mly heam“^, respectively. Dashed contours are for negative levels. The 
maps are corrected for the primary heam response. Shown in the bottom right-hand corners are the synthesized beams: FWHM = (/.'SI x Cf.'46 with 
RA. = -11°, FWHM = 0'.'33 x 0'.'29 with RA. = -21°, and FWHM = 0'.'33 x 0'.'29 with RA. = -60°, respectively. Rlus sign: 3.6 cm continuum 
source VLA 9 (Figure 3). Filled circle: CH 3 OH 6.7 GHz maser source (Green et al. 2012). The size of each marker represents the beam size of the 
corresponding observations. 


The 3.6 cm continuum image from the AR 465 observations 
shows that VLA 9 has an extended structure (Figure 3). The 
peak position is 02000 = 05h41“44!137, ^2000 = -01°54'46"06, 
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Fig. 3.— Map of the 3.6 cm continuum emission toward the FIR 
4 region, from the AR 465 data set. The contour levels are 1, 2, 4, 
and 8 x 0.03 mjy beam~^, and the rms noise is 0.011 mjy beam“^. 
Dashed contours are for negative levels. The map is corrected for the 
primary beam response. Filled ellipse: synthesized beam. FWHM = 
0'.'25 X 0'.'23 with RA. = -23°. Open ellipse: effective beam (apparent 
width of a point source, taking the bandwidth smearing into account) 
at the position of VLA 9, determined from other unresolved sources in 
the map. FWHM = 0'.'33 x O'.'24 with RA. = -10°. Large plus sign: 
H 2 O maser spot 1 (Section 3.2). Orange filled circle: CH 3 OH 6.7 GHz 
maser source (Green et al. 2012). The size of each marker represents 
the beam size of the corresponding observations. The straight line in 
the bottom left-hand corner corresponds to 100 AU at a distance of 
415 pc. 


consistent with Rodriguez et al. (2003). The total flux density 
is 0.57 ± 0.02 mJy. The source is elongated in the northwest- 
southeast direction. The extent of elongation is larger than 
the effective beam (synthesized beam degraded owing to the 
bandwidth smearing). In addition, the intensity distribution is 
slightly curved and asymmetric with respect to the peak po¬ 
sition. An elliptical Gaussian fit gives a deconvolved size of 
FWHM = 0'.'37 x (y.'09 with a position angle (RA.) = -25° (de¬ 
convolved with the effective beam). The major axis corresponds 
to FWHM= 150 AU. 

The spectrum of VLA 9 in the 6.2-3.6 cm interval is flat (Ta¬ 
ble 2), which suggests that the centimeter continuum is optically 
thin free-free emission. This spectrum rules out the possibility 
of VLA 9 being an ultracompact/hypercompact HII region. The 
elongated structure and the flat spectrum suggest that VLA 9 is 
a radio thermal jet. The 2.0 cm flux density is higher than those 
of the longer wavelengths, and the spectral index in the 3.6-2.0 
cm interval is ^1.2 (Figure 4), which suggests that a fraction of 
the 2.0 cm flux may come from dust. 

3.2. Water Maser 

The image from the AF 354 data set revealed detectable H 2 O 
maser spots near FIR 4 and FIR 6n. The FIR 6n maser source 
was already reported by Furuya et al. (2003), and we will omit 
descriptions of this source. (Also see the clarification of source 
identification in Choi et al. 2012a.) No maser source was de¬ 
tected from the AM 749 data set. The 3tT detection limits are 8 
mJy beam“^ at FIR 6n and 19 mJy beam“^ at FIR 4. The AM 
780 data set revealed detectable H 2 O maser spots near FIR 6n 
only. The detection limit at FIR 4 is 40 mJy beam“^. The results 
are summarized in Table 3. 

Three H 2 O maser spots associated with FIR 4 were detected 
from the AF 354 observations: two weak sources (spots 1 and 
2) near the systemic velocity of the ambient cloud (Vlsr = 1 LO 
km s“^; Schulz et al. 1991) and a strong one (spot 3) at a red- 
shifted velocity (Figures 5 and 6). The velocity difference be¬ 
tween them is ^8 km s“^. The H 2 O maser source detected by 
the single-dish observations in 2012 had a velocity of Flsr = 
13 km s“^ (Choi et al. 2012a). 
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Table 3 

NGC 2024 FIR 4/6 WATER MASER SOURCE PARAMETERS 


Region 

Epoch 

Spot 

Peak Position^ 

Q:j2000.0 5 j2000.0 

Centroid Velocity 
(km s“^) 

Peak Flux Density ° 
(Jy beam”^) 

FIR 4 

1999 Feb 26 

1 

05 41 44.12 

-01 54 46.4 

11.3 

0.16 ±0.04 



2 

05 41 44.18 

-01 54 48.9 

11.4 

0.10 ±0.04 



3 

05 41 44.02 

-01 54 43.3 

19.1 

4.29 ±0.04 

FIR 6n 

2003 Sep 15 

4 

05 41 45.155 

-01 56 00.64 

10.7 

3.536 ± 0.009 



5 

05 41 45.171 

-01 56 00.60 

13.7 

0.897 ± 0.009 


17.1 1.546 ±0.009 

21.9 0.987 ± 0.009 


® Units of R.A. are hours, minutes, and seconds, and units of decl. are degrees, arcminutes, and arcseconds. 
^ Flux density at the peak channel, corrected for the primary beam response. 



og u (GHz) 


Fig. 4.— Spectral energy distributions of the sources detected in the 
three wavelength bands (Table 2). The flux uncertainties (including the 
uncertainty of the flux scale) are smaller than the sizes of the markers. 


A Gaussian fit to spot 3 shows that it is slightly larger than a 
point source, with a deconvolved size of FWHM = I'.'O x 0'.'2. 
The intensity distribution of spot 3 can be described as a com¬ 
bination of a strong point-source component and a weak (below 
^10% of the peak) extended component. The weak component 
is elongated in the north-south direction (P.A. = -10°), which 
caused the peak at ^^19 km s“^ in the spectrum of spot 1 (Fig¬ 
ure 6(a)). Since this low-level elongation of spot 3 is symmetric 
around the intensity peak and has a position angle exactly the 
same as that of spot 3 with respect to the phase-tracking cen¬ 
ter, it is probably an artifact of the imaging process, not a real 
source structure. It is probably owing to the position far from 
the phase-tracking center and the velocity near the edge of the 
bandpass. 



RA (J3000) 


Fig. 5.— Maps of the H 2 O maser emission in the FIR 4 region. The 
maser spot numbers are labeled. The violet contours show the intensity 
integrated over Ulsr = (10.7, 11.7) km s“^, and the contour levels 
are 0.5 and 0.75 times the peak intensity of spot 1 (0.14 Jy beam“^ 
km s“^). The red contours show the intensity integrated over (18.6, 
19.6) km s“^, and the contour levels are 0.25, 0.5, and 0.75 times the 
peak intensity (2.67 Jy beam“'^ km s“^). The maps are corrected for 
the primary beam response. Shown in the bottom right-hand comer is 
the synthesized beam: FWHM = 3'.' 1 x 2'.' 1 with a P.A. = -84°. Gray 
dashed contours: map of the 4.5 pm emission from the Spitzer data 
archive. The contour levels are 1, 2, and 4 x 200 MJy sr~^. Plus sign: 
3.6 cm continuum source VLA 9 (Figure 3). Red arrow: direction of 
the redshifted CO outflow (Chandler & Carlstrom 1996). 


4. DISCUSSION 

The positions of the radio continuum source and the I-I 2 O 
maser spots are shown in Figure 7, together with other com¬ 
pact objects in the FIR 4 region. Most of them are distributed 
along a straight line in the northwest-southeast direction, which 
suggests that they belong to a single protostar-outflow system. 
The extended infrared feature in the 4.5 pm map (Figure 5), ex¬ 
tending from the infrared peak position toward the northwest, 
can also be seen in the Spitzer maps of other bands (3.6, 5.8, 
and 8.0 pm), suggesting that it is seen in continuum emission. 
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Fig. 6.— Spectra of the H 2 O maser line in the FIR 4 region, (a) Spec¬ 
trum of spot 1. The peak at ~19 km s“^ is owing to the emission from 
spot 3. (b) Spectrum of spot 2. (c) Spectrum of spot 3. Each spec¬ 
trum was taken at a single pixel at the peak position of the correspond¬ 
ing maser spot. Vertical dotted line: velocity of the ambient dense gas 
(Vlsr = 11.0 km s"^ Schulz et al. 1991). 



Fig. 7.— Compact objects in the FIR 4 region. Small plus sign: 3.6 
cm continuum source VLA 9 (the peak position in Figure 3). Large 
plus signs: H 2 O maser spots (Figure 5). Double circle: near-IR source 
(2MASS Point Source Catalogue). Open circles: peak positions of the 
Spitzer maps (green: 4.5 ^m, red: 8.0 ^m). Gray square: 3 mm contin¬ 
uum source (Eisner & Carpenter 2003). Crosses: X-ray sources (Skin¬ 
ner et al. 2003). The size of each marker represents the beam size of 
the corresponding observations for radio sources or the pixel size for 
infrared and X-ray sources. 


The position and morphology of this feature agree with those of 
the reflection nebulosity seen in the H and K bands (Moore & 
Chandler 1989). 


4.1. FIR 4 Protostar 

The nature of the YSO embedded in the FIR 4 cloud core 
has been controversial (Mezger et al. 1988, 1992; Moore & 
Chandler 1989; Minier et al. 2003; also see Section 4.2 of Choi 
et al. 2012a). Close examinations of the available information 
reveal that a major part of the confusion stems from the in¬ 
frared/millimeter source position issue. 

Moore & Chandler (1989) considered that the point-like 
near-infrared (near-IR) source represents the embedded YSO. 
However, peak positions of the near-IR to mid-infrared (mid- 
IR) maps do not coincide with each other (Figure 7). They 
move toward the southeast with increasing wavelength, along 
the outflow axis. They tightly follow the straight line defined by 
VLA 9 and the 2MASS source position. This position shift sug¬ 
gests that the detected near-IR emission comes from the scat¬ 
tered light escaping through the outflow cavity. The near-IR po- 
larimetry is consistent with this interpretation (Kandori et al. 

2007) . This finding solves the extinction problem pointed out 
by Mezger et al. (1992). They argued that the expected extinc¬ 
tion of the FIR 4 cloud core is too high for any embedded YSO 
to be directly detectable in near-IR. The central protostar is in¬ 
deed invisible in near-IR. Judging from the trend of the infrared 
position shift, the true position of the protostar may lie some¬ 
where to the southeast of the 8.0 /rm peak position. The best 
candidate is the peak position of the radio thermal jet VLA 9. 

High-resolution interferometric imaging in the millimeter 
continuum can provide additional information on the position 
of the protostar. Known from interferometric observations are 
the 3 mm source positions reported by Chandler & Carlstrom 
(1996) and Eisner & Carpenter (2003). Their positions of FIR 
4 are inconsistent and different by ~2'.'2 (comparable to their 
beam sizes). A comparison of the positions of FIR 5 and 6 in 
the literature (Wiesemeyer et al. 1997; Lai et al. 2002; Alves et 
al. 2011; Choi et al. 2012b) shows that the positions of Eisner & 
Carpenter (2003) are accurate to within the beam size, but those 
of Chandler & Carlstrom (1996) are systematically shifted with 
respect to those of the other references. The amount of shift 
is ^3" toward the northwest. This astrometric problem caused 
some confusion in source identifications (Minier et al. 2003; 
Skinner et al. 2003). The EIR 4 position of Eisner & Carpen¬ 
ter (2003) agrees with VLA 9 within the beam size (Eigure 7). 
Therefore, VLA 9 is the best indicator of the position of the EIR 
4 protostar, and there is no evidence for multiple YSOs in the 
EIR 4 dense core. 

Setting aside the CH 3 OH maser issue (Section 4.2), most of 
the evidence indicates that EIR 4 contains a low-mass proto¬ 
star. Mezger et al. (1992) derived 10 Mq for the mass of the 
cloud core assuming that there is no internal heating and the 
core is cold (~19 K). However, this assumption is inconsistent 
with recent findings. Later studies showed that the temperature 
is higher (50-80 K; Mangum et al. 1999; Watanabe & Mitchell 

2008) . Then the mass of the EIR 4 core may be much smaller: 
0.8-1.7 Mq (Chandler & Carlstrom 1996; Visser et al. 1998; 
Watanabe & Mitchell 2008). Considering that EIR 4 is in an 
early stage of protostellar evolution, the core mass seems to be 
too small to form a high-mass/intermediate-mass star. 

Moore & Yamashita (1995) estimated a bolometric luminos¬ 
ity of ~25 Lq for the YSO in EIR 4, but this value is uncertain 
because the far-infrared (far-IR) flux was not included in the 
calculation. Lis et al. (1991) estimated a luminosity of 3-7 Lq 
for the internal heating source (or much smaller if the external 
radiation field is highly enhanced). Therefore, the luminosity 
also indicates that the EIR 4 core contains a low-mass protostar. 
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One of the major difficulties in understanding the nature of 
the FIR 4 protostar is that the nebulosity associated with the 
ionization front and the infrared-bright stars in the nearby re¬ 
gion prevents accurate photometry in the far-IR band where the 
spectral energy distribution peaks. The luminosity from avail¬ 
able data (discussed in the previous paragraph) and the lack of 
an ultracompact/hypercompact H II region firmly rule out the 
likelihood of FIR 4 containing a high-mass YSO. The FIR 4 
protostar is most likely a low-mass object, but the possibility of 
an intermediate-mass (3^ Mq) object cannot be clearly ruled 
out. Further investigations are needed to settle this issue. 

4.2. Methanol Maser Source 

One of the interesting features of FIR 4 is the CH 3 OH 6.7 
GHz maser emission (Minier et al. 2003). This class II maser 
is considered to be exclusively associated with massive star¬ 
forming regions (Menten 1991; Xu et al. 2008). The CH 3 OH 
maser spot of FIR 4 is one of the nearest class II maser sources 
known (Green et al. 2012) and provides an excellent opportunity 
to investigate the physical relation between YSOs and maser 
sources. 

The CH 3 OH maser source is closely associated with VLA 9 
and H 2 O maser spot 1 (Figure 3). The position marked in Figure 

3 was taken from Figure 4 of Green et al. (2012), which is more 
accurate than the one given by Minier et al. (2003). The velocity 
of the CH 3 OH maser line (V):,SR = 12.3 km s“^; Minier et al. 
2003; Green et al. 2012) is redshifted by 1.3 km s“^ relative to 
the systemic velocity (11.0 km s“^; Schulz et al. 1991) or by 1.0 
km s“^ relative to the velocity of the H 2 O maser spot 1 (Table 
3). The widths of the CS and H 2 CO lines from the FIR 4 dense 
core are ~2 km s“^ (Schulz et al. 1991; Mangum et al. 1999). 
The line profiles of H 2 CO are nearly symmetric and have no 
line wing (Mangum et al. 1999), which suggests that the line 
width mostly reflects the turbulent motion of dense molecular 
gas. Therefore, it is unclear if the velocity shift of the CH 3 OH 
maser is related to the outflow, disk rotation, or any other motion 
in the cloud core. 

The angular separation between the CH 3 OH maser spot and 
the peak of the VLA 9 is ~(y.'21. (The beam size of the CH 3 OH 
line observations of Green et al. (2012) is ^O'.'OT. The position 
accuracy of the VLA is usually better than O'.' 1.) The separation 
may not be completely accurate because the source positions 
came from different observations. Simultaneous observations in 
both the continuum and the CH 3 OH maser line are necessary to 
measure the separation more accurately. 

Minier et al. (2003) argued that FIR 4 contains a high-mass 
or intermediate-mass protostar, based on the assumption that the 
CH 3 OH maser action may require warm (> 175 K) gas heated 
by the protostar. They thought that the CHsOH-rich gas in a 
low-mass star-forming region exists far (^ 100 AU) from the 
central YSO where the temperature is too low to drive a maser 
action. The proximity of the CH 3 OH maser to VLA 9 suggests 
that this argument needs to be reexamined. 

The usual assumption is that the maser spot is located in 
the heated envelope of a YSO (Minier et al. 2003). If the FIR 

4 protostar is located exactly at the peak position of VLA 9, 
the projected separation between the CH 3 OH maser source and 
the protostar is ~90 AU. The dust temperature at the CH 3 OH 
maser-emitting region can be estimated using Equation (2) of 
Motte & Andre (2001). Assuming that the distance from the 
protostar to the maser spot is ~120 AU (considering the pro¬ 
jection effect) and that the protostellar luminosity is 3-25 Lq 
(Section 4.1), the dust temperature would be 40-70 K. How¬ 
ever, models of class II masers suggest that mid-IR photons 


from warm (100-200 K) dust are necessary for radiative pump¬ 
ing (De Buizer et al. 2000; Gragg et al. 2005). Then the CH 3 OH 
maser of FIR 4 suggests several possibilities. (1) The class II 
maser action may operate at temperatures lower than the theo¬ 
retically expected range. (2) The FIR 4 maser spot may have an 
extra source of heating, such as the outflow or the external radi¬ 
ation field. (3) The maser spot may be on the circumstellar disk 
irradiated by photons coming directly from the protostar. High- 
resolution (~ 0'.'2 or better) observations in the dust continuum 
may be helpful in understanding the relation between the pro¬ 
tostar and the CH 3 OH maser source. The flat radio spectrum, 
small mass, and low luminosity (Section 4.1) rule out the exis¬ 
tence of a high-mass YSO in the FIR 4 region, and the CH 3 OH 
maser must be excited by the FIR 4 low-mass protostar regard¬ 
less of the detailed heating and pumping mechanisms. 

4.3. FIR 4 Outflow 

The three H 2 O maser spots are distributed along a straight 
line in the northwest-southeast direction with a RA. of about 
-23° (Figure 5). This line passes through VLA 9 and agrees 
with the axis of the outflow in this region: the northwestern 
lobe indicated by the reflection nebulosity and the southeastern 
lobe traced by the CO line (Moore & Chandler 1989; Chan¬ 
dler & Carlstrom 1996). The unipolar CO outflow is redshifted 
(Chandler & Carlstrom 1996). The southeastern CO outflow 
may be on the far side of the FIR 4 core, and the northwestern 
infrared nebulosity (outflow cavity) may be on the near side. In¬ 
terestingly, the H 2 O maser spot 3 in the northwest is redshifted, 
which may reflect the velocity of the local shocked gas. 

The centimeter continuum source VLA 9 (Figure 3) is lo¬ 
cated close to H 2 O maser spot 1. VLA 9 is elongated, and its 
major axis agrees with the direction of the outflow, which sug¬ 
gests that VLA 9 is a radio thermal jet at the base of the bipolar 
outflow. With respect to the intensity peak, the emission struc¬ 
ture is more extended toward the northwest than toward the op¬ 
posite direction. 

One of the X-ray sources reported by Skinner et al. (2003) 
is located close to H 2 O maser spot 3 (Figure 7). The separation 
is ~0'.'7 or 300 AU. The H 2 O maser and the X-ray emission 
are not necessarily coming from the same volume of gas, but 
the northwestern outflow of FIR 4 may be responsible for both 
phenomena. Only a few young protostars are known to show 
the H 2 O maser and X-ray emission together. The YSOs in the 
NGC 2071 cluster (IRS 1/3 and VLA 1) are notable examples 
(Torrelles et al. 1998; Skinner et al. 2009; Trinidad et al. 2009). 
High-resolution observations in the H 2 O maser line and deeper 
observations in the X-ray are necessary to understand the phys¬ 
ical relation between the maser-producing shocked gas and the 
X-ray-emitting hot plasma. 

The four X-ray photons detected in the FIR 4 region were 
hard (~5.9 keV; Skinner et al. 2003), and the emission mecha¬ 
nism is difficult to explain. Skinner et al. (2003) suggested that 
the emission may be caused by stellar flares, and additional pos¬ 
sibilities (mainly relevant to NGC 2071 IRS 1) are listed in Sec¬ 
tion 5.2 of Skinner et al. (2009). The large separation (^2'.'8) be¬ 
tween the X-ray source and the protostar (VLA 9) seems to rule 
out all of these possibilities but the outflow model. However, 
if the X-ray-emitting plasma is produced by the outflow shock¬ 
ing on dense ambient gas, the required shock speed is ~2200 
km s“^ (see the discussion in Section 5.2.1 of Skinner et al. 
2009). This speed is unreasonably high for a YSO outflow. Yet 
another model was proposed by Lopez-Santiago et al. (2013) to 
explain the X-ray emission of HH 80. In this model the hard 
X-ray emission is synchrotron radiation produced by nonther- 
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Fig. 8.— Maps of the 3.6 cm continuum emission around VLA 1 and VLA 21, from the AR 465 data set. The contour levels are 1, 2, 4, 8, 16, 
32, and 64 X 0.05 mjy beam The maps are corrected for the primary beam response. Shown in the bottom right-hand corner is the synthesized 
beam: FWHM = O'.' 25 X o'.' 23 with P.A. = —23°. The elongated shape of VLA 1 is owing to the bandwidth smearing. Dashed contours: maps of the 
4.5 /rm emission from the Spitzer data archive. The contour levels are 1, 2, and 4 x 200 MJy sr“^ for (a) and 1, 2, 4, 8, 16, and 32 x 20 MJy sr~^ 
for (b). Cross: X-ray source SGB 193 (Skinner et al. 2003). 


mal processes in the magnetized jet. It is unclear if this model 
applies to the case of FIR 4, because there is no detection of 
synchrotron radiation or strong magnetic fields. Therefore, the 
origin of the hard X-ray remains a puzzle. 

5. NEARBY OBJECTS 

5.1. Continuum Sources 

Rodriguez et al. (2003) detected a cluster of 25 radio sources 
in the NGC 2024 region. Based on a statistical consideration, 
they suggested that practically all of these sources are physi¬ 
cally associated with NGC 2024. Many of these radio sources 
are indeed associated with infrared sources and other star for¬ 
mation activities such as outflows. The nature of each source 
can be examined with the spectrum of the centimeter contin¬ 
uum from ionized gas (Reynolds 1986; Anglada et al. 1998). 
Optically thin free-free emission (e.g., from a thermal radio jet) 
shows a small spectral index (ct < 1). Optically thick free- 
free emission (e.g., from a high-density H II region) shows a 
larger spectral index. (At shorter wavelengths, the emission can 
become optically thin, and the spectrum can become flat.) Non- 
thermal emission (e.g., from a magnetic flare) usually (but not 
always) shows a negative spectral index (Giidel 2002). Table 2 
lists the spectral indices of 10 radio sources that were detected 
in both 6.2 and 3.6 cm maps. They can be classified into three 
categories depending on the spectral slope. The type of radio 
source is listed in the last column of Table 2. (Also see the de¬ 
scriptions of several sources in Rodriguez et al. 2003.) 

5.1.1. Positive-spectrum Sources 

Four sources (VLA 16, VLA 19, VLA 21, and VLA 25) have 
relatively large spectral indices. The centimeter continuua of 
these sources (except VLA 21) are probably (partially) optically 
thick free-free emission. 

VLA 16 is elongated in the northwest-southeast direction 
(see Figure 5 of Rodriguez et al. 2003). An elliptical Gaussian 
fit gives a deconvolved size of FWHM = 0'.'38 x 0'.'15 with 
RA. = -39°. The nature of VLA 16 is unclear. It has an X-ray 


counterpart but remains undetected in the infrared (Skinner et 
al. 2003). Rodriguez et al. (2003) suggested that VLA 16 forms 
a binary system with VLA 15 (IRS 2b). 

VLA 19 (IRS 2) is the brightest radio/infrared source in this 
region. The flux density increases steadily with the frequency, 
and the spectral index in the 6.2-2.0 cm interval is ~ 1.0 (Figure 
4). Kurtz et al. (1994) listed VLA 19 (G206.543-16.347) as an 
ultracompact H II region containing a B2 star. Lenorzer et al. 
(2004) suggested that the radio emission comes from a stellar 
wind of an early-B-type star, recombining at a radius of ~100 
AU. 

VLA 21 is associated with an infrared source and an X- 
ray source (Figure 8). Based on the high flux variability in 
a timescale of a few days, Rodriguez et al. (2003) suggested 
that the 3.6 cm continuum is gyrosynchrotron emission, which 
seems to conflict with the positive spectral index. (Considering 
the flux variability, however, the uncertainty of the spectral in¬ 
dex can be much larger than the value given in Table 2.) A pos¬ 
sible explanation may be that the radio flux of VLA 21 has both 
thermal and nonthermal components. The quality of the multi- 
frequency data is not good enough to understand this problem. 
Sensitive monitoring observations are necessary to understand 
the variability of the flux and the spectral index. 

5.1.2. Flat-spectrum Sources 

Four sources (VLA 1, VLA 3, VLA 9, and VLA 24) have 
relatively small spectral indices. The centimeter continuum of 
these sources (except VLA 24) are probably optically thin free- 
free emission. 

Among the eight YSOs listed in Table 2, only VLA 9 (FIR 
4) is associated with a dense molecular cloud core (Chandler & 
Carlstrom 1996; Eisner & Carpenter 2003). VLA 9 is probably 
the youngest object among these YSOs. 

Based on the circular polarization and mild flux variability 
Rodriguez et al. (2003) suggested that VLA 24 may be a young 
low-mass star showing gyrosynchrotron emission. The positive 
spectral index suggests that the radio flux of VLA 24 may also 
be a mixture of thermal and nonthermal components. 





5.1.3. Negative-spectrum Sources 

Two sources (VLA 2 and VLA 8) have significantly nega¬ 
tive spectral indices. The centimeter continuum of these sources 
may be nonthermal emission. 

VLA 8 is elongated in the north-south direction (see Figure 3 
of Rodriguez et al. 2003). An elliptical Gaussian fit gives a de¬ 
convolved size of FWHM = (y.'68 x 0'.'31 with a P.A. = -6°. The 
spectral index in the 6.2-2.0 cm interval is about -1.5 (Figure 
4). Based on the morphology and the orientation of the elon¬ 
gation Rodriguez et al. (2003) suggested that VLA 8 may be 
an ionized proplyd. However, the radio proplyd model conflicts 
with the steeply decreasing radio spectrum. Considering the ex¬ 
tended structure, no day-scale variability, and lack of infrared 
counterpart, it is unlikely that VLA 8 is a magnetically active 
YSO showing gyrosynchrotron emission. VLA 8 is probably a 
background object such as an extragalactic radio jet. 

5.2. FIR 6n 

Two H 2 O maser spots associated with FIR 6n were detected 
from the AM 780 observations: one near the systemic velocity 
of the ambient cloud (Vlsr = 10.8 km s“^; Choi et al. 2012a) 
and the other at redshifted velocities (Figure 9). The redshifted 
source (spot 5) shows three velocity components (Figure 10). 
The H 2 O maser sources detected previously (Furuya et al. 2003; 
Choi et al. 2012b) positionally coincides with spots 4/5, but 
their beam sizes were too large to separate the two spots. 
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Fig. 9.— Maps of the H 2 O maser and 3.6 cm continuum emission 
in the FIR 6n region. The violet contours show the H 2 O intensity in¬ 
tegrated over Vlsr = (9.3, 12.0) km s“^, and the contour levels are 
0.25, 0.5, and 0.75 times the peak intensity (8.88 Jy beam“^ km s“^). 
The red contours show the H 2 O intensity integrated over (12.0, 23.8) 
km s“^, and the contour levels are 0.25, 0.5, and 0.75 times the peak 
intensity (7.13 Jy beam“^ km s“^). The maser spot numbers are la¬ 
beled. The gray contours show the 3.6 cm continuum intensity from 
the AR 465 data set, and the contour levels are 1, 2, 4, and 8 x 0.03 
mjy beam“^. The maps are corrected for the primary beam response. 
Filled ellipse: synthesized beam of the H 2 O maps. FWHM = C/.'l3 x 
OiHl with RA. = -17°. Open ellipse: synthesized beam of the 3.6 cm 
continuum map. FWHM = O'/25 x 0'.'23 with P.A. = -23°. The straight 
line in the bottom left-hand corner corresponds to 100 AU at a distance 
of 415 pc. 
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Fig. 10.— Spectra of the H 2 O maser line in the FIR 6n region, (a) 
Spectrum of spot 4. (b) Spectrum of spot 5. Vertical dotted line: veloc¬ 
ity of the ambient dense gas (Vlsr = 10.8 km s”*^; Choi et al. 2012a). 


Figure 11 shows the positions of the H 2 O maser spots and 
other objects in the FIR 6 region. The maser spots coincide 
with the millimeter continuum source FIR fin. The separation 
between the two spots is much smaller than the angular reso¬ 
lutions of the millimeter continuum maps available (Lai et al. 
2002; Alves et al. 2011; Choi et al. 2012b). The centimeter con¬ 
tinuum source VLA 14 (Rodriguez et al. 2003) coincides with 
H 2 O maser spot 5 (Figure 9). VLA 14 is unresolved. The peak 
position and flux density given by Rodriguez et al. (2003) are 
consistent with those from our analysis. 

The two maser spots are separated in the east-west direction 
(Figure 9). The separation (0'/24) corresponds to 100 AU, and 
the relative position angle is 81°. This value is similar to the 
position angle of the 30" scale bipolar outflow traced by the CS 
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Fig. 11.— Positions of the compact objects in the FIR 6 region. Filled 
circles: H 2 O maser spots (Figure 9). Grayscale: map of the 6.9 mm 
continuum emission (see Figure 3 of Choi et al. 2012b). Dashed con¬ 
tours: map of the 4.5 pm emission from the Spitzer data archive. The 
contour levels are 2, 3, 4, 5, 6, 7, and 8 x 13 MJy sr“^. The emission 
in the northwestern comer of the map comes from an infrared star. 
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line and the high-velocity CO J = 1 0 emission (~75°; Chan¬ 

dler & Carlstrom 1996). The extended infrared feature in the 4.5 
fim map (Figure 11) shows two emission peaks, one in the east 
and the other in the west, and seems to be related to the out¬ 
flow. The 4.5 /xm emission feature is not obvious in the Spitzer 
maps of the other bands, suggesting that it probably comes from 
line emission in the 4.5 pm band. The separation between the 
4.5 pm peaks is ~10", and the relative position angle is ~106°, 
which is different from that of the CS/CO outflow by ~30°. 
Alves et al. (2011) showed that the outflow traced by the CO J 
= 3—^2 line has a complicated morphology and suggested that 
the outflows of FIR 6 n and FIR 5 may be interacting. 

FIR 6 n is the most active H 2 O maser source in the whole 
NGC 2024 molecular ridge region (Choi et al. 2012a). The na¬ 
ture of the YSO in FIR 6 n is unclear because, in the millime¬ 
ter/infrared wavelengths, FIR 6 n is often either undetectable or 
affected by FIR 6 c. Considering that FIR 6 n drives a strong 
outflow and is invisible in the near-IR, the FIR 6 n cloud core 
may contain a low-mass protostar. The nature of the stronger 
millimeter/submillimeter source FIR 6 c is even more unclear. 
There is no molecular outflow associated with FIR 6 c, but a 
redshifted H 2 O maser was detected (Choi et al. 2012b). FIR 6 c 
was undetected in the 3.6 cm continuum. The spectral energy 
distribution of FIR 6 c rises steeply around 6.9 mm and is flat 
around 1 mm, which suggests that FIR 6 c probably contains a 
hypercompact H II region (Choi et al. 2012b). If so, the central 
object may be an early-B type star. High-resolution imaging in 
the millimeter continuum is necessary to understand this pecu¬ 
liar object. 

6 . SUMMARY 

The star formation activities in the NGC 2024 region were 
investigated by analyzing several archival data sets from the 
VLA. The nature of centimeter continuum sources was exam¬ 
ined using a data set of multifrequency observations. The source 
positions and structures were inspected using the data set of 3.6 
cm continuum observations previously published by Rodriguez 
et al. (2003). Outflow activities were studied using several data 
sets of H 2 O maser observations. Careful examinations of the 
positions of compact objects in the FIR 4 region provide a sim¬ 
ple and comprehensible picture of the star formation activities 
around the FIR 4 protostar. The main results are summarized as 
follows. 

1. The centimeter continuum source VLA 9, the source as¬ 
sociated with FIR 4, shows a flat spectrum in the 6 .2-3.6 cm 
interval. VLA 9 is elongated in the direction consistent with the 
FIR 4 outflow. VLA 9 is a radio thermal jet and the best indica¬ 
tor of the protostellar position. The flat spectrum rules out the 
possibility of a high-mass YSO in the FIR 4 region. Previously 
known mass and luminosity estimates suggest that the FIR 4 
protostar is most likely a low-mass object. 

2. The CH 3 OH 6.7 GHz maser source in FIR 4 is located 
close to VLA 9. The separation between them is ~(y.'2. They 
must have a close physical relation, but the detailed amplifica¬ 
tion mechanism of the maser is unclear. The CH 3 OH maser of 
FIR 4 demonstrates that the class II maser phenomenon does 
not necessarily require the presence of a high-mass YSO. 

3. The H 2 O maser spots of FIR 4 are distributed along the 
FIR 4 outflow axis. The strongest spot is associated with an X- 
ray source. The FIR 4 outflow is a rare example of protostellar 
outflows showing the H 2 O maser and X-ray emission together. 
Further studies are needed to understand the origin of the hard 
X-ray emission. 


4. The properties of other radio sources in the field of view 
were investigated. VLA 1, VLA 3, and VLA 25 (IRS 4) may 
be radio thermal jets driven by YSOs. VLA 2, VLA 21, and 
VLA 24 may be magnetically active YSOs. VLA 8 may be a 
background object emitting nonthermal radiation. VLA 16 is a 
radio source of unknown nature. VLA 19 (IRS 2) may be an 
early-B-type star. 

5. The two H 2 O maser spots of VLA 14 (FIR 6n) are dis¬ 
tributed along an east-west line that is consistent with the axis 
of the FIR 6n molecular outflow. The FIR 6n cloud core may 
contain a low-mass protostar. 
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